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Abstract 

In this article we exploit the known commutative family in Y(g[ n ) - the Bethe 
subalgebra - and its special limit to construct quantization of the Gaudin inte- 
grable system. We give explicit expressions for quantum hamiltonians QIk(u), k = 
1, . . . , n. At small order k = 1, . . . , 3 they coincide with the quasiclassic ones, even 
in the case k = 4 we obtain quantum correction. 
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1 Introduction 



Quantization. The quantization problem has different aspects and different approaches 
therein. The quantization consists in constructing a deformation of a Poisson algebra of 
functions on a phase manifold and an appropriate representation. We are interested 
in the first part of quantization problem in a specific case of integrable systems. An 
integrable system is a hamiltonian system on a Poisson manifold (M 2n , it) where it is a 
Poisson bivector, given by a hamiltonian H such that there exists a family of algebraically 
independent functions Ii, . . . ,I n in involution {ij, Ij} = and H — l x . In this language 
an integrable system is a subalgebra 3 C C(M) which is also commutative subject to the 
Poisson bracket and is generated by the appropriate number of generators. This property 
of being appropriate is even more ambiguous at the quantum invariant level, so we do not 
discuss it here. 

Traditionally the quantization of a Poisson algebra on C(M) is an associative algebra 
A(M) isomorphic as a linear space to the space of formal series C(M)[h] with a fixed 
isomorphism of linear spaces 

q : C(M)[H] -> A(M) 
called a quantization correspondence, such that 

q(f)*q(9)=q(f9)+0(h) and q(f) * q(g) - q(g) * q(f) = %({/, g}) + 0(h 2 ) 

where /, g are functions on M, * is the multiplication on A(M) and {/, g} is the Poisson 
bracket. The quantization of an integrable system in this context is such a quantization 
of a Poisson algebra that q($s) is a commutative subalgebra (not only at the first order 
in h) in A(M). The importance of this condition is obvious, as like as in the classical case 
the dynamics can be restricted to the common level of integrals, the quantum problem of 
finding the eigenvectors can be restricted to the common eigenspace of quantum integrals. 

The main example of this paper is the quantization of the space g* for a Lie algebra g = 
gl n with the Kirillov-Kostant Poisson bracket. The space of functions in this case is just 
the symmetric algebra C(M) = S*(g). The quantum algebra is the universal enveloping 
algebra U (g) deformed in a standard manner: one takes new generators ejj = he^ than 
the defining relations in U(g)[h] became 

[eij,e k i] = h(S jk eu - 8ue kj ). 

The quantization morphism on homogeneous polynomials is given by fixing some normal 
ordering. The classical limit consists in tending h \— > 0. 

The same procedure can be realized without thinking of formal parameter due to the 
structure of filtered algebra on U(gl n ). Indeed, the quantization can be considered as a 
map from S*(gi n ) to just U(gl n ) given by the same normal ordering. Let us note that 
the choice of such an ordering does not affect the result of passing to the associated 
graded algebra Gr(U(gl n )) = S*(gi n ). Hence one can consider the factorization morphism 
p : U(gl n ) — > Gr(U(gl n )) as a classical limit. Indeed, the composition of quantization and 
factorization is an identity on homogeneous element of the symmetric algebra S*(gl n ). 
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Gaudin system. Let us consider the direct sum of k coadjoint jj[ n -orbits @iOi with 
an induced Poisson (symplectic) structure and the diagonal symplectic GL n -action. An 
element of this space can be represented by the Lax operator 

£ <"> = £^ (i) 

i 

with <i>j e Oi. Fixing a basis one can think of <£>j as a matrix valued function on Oi with 
matrix elements fki - generators of gl n considered as linear functions on the i-th orbit. 
Then the functions 

Ii(u) = TrU{u) commute {Ii(u), Ij(v )} = 

and are invariant under the diagonal GL n -action. The images of these functions on the 
symplectic reduction space M = @iOi/ /GL n give the integrable system. The proof of the 
integrability can be found for example in PQ and is based on the r-matrix representation 
of the Poisson bracket: 

P 

{L(u) ® L(v)\ = [r(u - v),L(u) ® 1 + 1 ® L(v)\ with r(u) = — 

u 

where P is a permutation operator in the tensor product. There exist several partial 
responses to the quantization problem of this integrable system. For example in the gl 2 
case just the same expressions 

h{u) = TrL 2 {u) 

with fki changed to the U(gl z ) generators give the commutative family of elements in 
U(gl z )® k . The attempt to generalize this result to gl n in [2] was only partially successful. 
There was obtained that quantum expressions for h(u), h(u) commute but it is not true 
in general for the higher integrals. 



Main result. Consider the quantum Lax operator given by the expression Q where 
$j are matrices with elements in U(gl n )® k such that is the ej~i generator of gi n lying 

in the i-component of the tensor product. Consider the following differential operators 
applied to the function 1: 

Qh{u) = Tr^A^L^u) - d u )(L 2 (u) - d u ) . . . (L t (u) - d u )l (2) 

where is the antisymmetrizer in Mat®™ n . 

Main Theorem. The quantities Qhiu) constitute a commutative family in U(gl n )® k (u) 
in the following sense 

[QIi(u),QIj(v)] = 0, Vi,j = l,...n, Vw,v; 

and their classical limit in the sense of associated graded algebra coincides with the family 
of classical Gaudin hamiltonians. 
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Plan of the paper. In section ETT1 we recall the construction of the Yangian Y(Ql n ) and 
the Bethe subalgebra in it, which is a maximal commutative subalgebra for general choice 
of parameters. Further, considering the k-th tensor power of the evaluation homomor- 
phism we obtain a commutative family in £7(f|l n ) . In section POl we extend the Bethe 
subalgebra to the commutative family of formal expressions in h. This formal family gives 
the standard Bethe subalgebra at h = 1, whereas we use its special limit at h = 0. We also 
slightly change a basis using only linear transformations. The main result of the paper is 
contained in section 1231 

Acknowledgement. The author would like to thank his colleagues in ITEP for numer- 
ous discussions and stimulating atmosphere, especially A. Chervov, L. Rybnikov, as well 
as S. M. Kharchev and S. M. Khoroshkin for useful remarks. 



2 Quantum Gaudin system 
2.1 Bethe subalgebra 

The Yangian was introduced by Drinfeld |3. and take an important role in describing 
rational solutions of the Yang-Baxter Equation. Y(gl n ) is generated by elements tffi 
subject to specific relations on the generating function T{u) with values in Y(gl r 
End(C n ) 

i,j k 

where E^ are matrices with 1 on the i, j-th place. The relations on this generating 
function involve the Yang i?-matrix 



v 



R( u ) = 1- -V£ 



u — ' 

hi 



and reads as the following equation 

R(z - u)T!(z)T 2 (u) = T 2 (u)T x (z)R{z - u) 
in End(C n ) m <g> Y(gl n )[z,u), where 

T 1 (z) = J2 E ij ® 1 ® Uj (z) , T 2 (u) = 1 ® E U ® ( u ) ■ 



1,3 »J 



There is a known realization of Y(Ql n ) in U(gl n ) i.e. a homomorphism p x : Y(gl n ) — > 

1 + -^)^®^ = / :l + -, (3) 



where e^- are generators of g[ n . The following expression also provide a realization of 
Yangian in U(gl n )® k for the given /c-tuple of complex numbers a = (zi,...,Zk) 

T» = T\u- z x )T 2 (u -z 2 )... T\u - z k ), 
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where T l (u — z{) is the realization given by (jHJ) with lying in the /-th component of 
U(Ql n )® k - Let p a be the corresponding mapping p a : Y(gl n ) — > U(gl n )® k . 

Lemma 1 p a is a homomorphism. This is equivalent to the relation 

R{z - u)T?{z)T%{u) = T?{u)T?{z)R{z - u) 
in End{C n 'f 2 ® U(gl n )® k [z,u}. 

Proof. Here we use the fact that T{(z) commutes with T k (u) if I ^ k. Indeed, 

R(z - u)Tl{z - Zl )T*(z -z 2 )... T k (z - z k )T\[u - z x )T 2 2 [u - z 2 ) . . . T 2 k (u - z k ) = 

R{z - u)Tl{z - Zl )Tl{u - Zl )T?(z -z 2 )... T x \z - z k )T%{u - z 2 ) . . . T k (u - z k ) = 

T 2 \u - Zl )Tl{z - Zl )R(z - u)T*{z -z 2 )... T k (z - z k )T*(u - z 2 ) . . . T 2 k (u - z k ) = 

et cetera... 
□ 

The simple way to obtain the commuting family in U(g[ n )® k is to take the image 
p a (Z(Y(g[ n ))) of the center of Yangian. The homomorphism p a is not surjective and that 
is why this commuting family is not trivial (contains not only constants). The center of 
Yangian is given by 

qdetT a (u) = sgn(a)tZ (1)jl (u) . . .t a a{n)>n {u - n+ 1) 

see for example Proposition 2.7 from jlj. Unfortunately this family is too small, T a (u) is 
a rational function on u with k simple poles, qdetT a (u) has kn simple poles, the corre- 
sponding residues play the role of formal quantum hamiltonians whereas the dimension 
of the classical phase space is of order kn 2 . There is another construction enlarging the 
family of central elements in Yangian up to a maximal commutative subalgebra (see 0,0 
or 2.14 of [7]): let C be an n x n matrix, T(u) the generating function for the elements 
of Y(gl n ), A n is the matrix of antisymmetrizer in C nlX,n and 

Tiiu) = ^ 1 ® • • • ® 1® E M ®1 ® • • • ® 1 ® tki(u) 

kl 

then the coefficients of 

Ti(u) = Tr A n Ti(u)T 2 (u - 1) . ..T^u - i + l)C i+1 ...C n i = l,...n 
generate the commuting family 

[n(u),Tj(v)] = 

which is maximal if the matrix C has simple spectra. In the realization T a (u) each Ti(u) 
is a rational function on u with ik simple poles, the total number of formal quantum 
conserved quantities is equal to k ^ n+ 2 ^ n which is approximatively half the dimension of 
the classical phase space. 
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2.2 Formal deformation 

Let us consider the canonical rescaling of the universal enveloping algebra obtained by 
multiplying structure constants by h. The corresponding gl n representation of Yangian 
T(u) looks as follows 

T(u) = 1 + — . 

u 

It solves YBE with the rescaled Yang R-matrix 

R(u) = 1 - —. 

u 

The construction of the commuting family is reproduced literally in this case: 



n{u) = (TrA n T?(u)T?(u - h) . . . T«(u - (i - l)h)C m ...C n ) i = l,...n. (4) 



We will use later the following obvious 
Lemma 1 Let two formal in h expressions 



M(h) = M ^ and N ( h ) = J2 Nihi 



with values in some associative algebra A commute, then their first coefficients also com- 
mute 

[M m , N n ] = 0. 

It is also true for the first non-central coefficients of formal series. 

Remark 1 Let us note that the classical Gaudin Lax operator can be recovered as a 
classical limit of T a : 

T a {u) = l + hL(«) +0(h 2 ) 

where 

is the first non-central coefficient of the operator-valued matrix. This expression lies in 
the first filtration component and after taking classical limit remains of the same form 
(with described in Introduction exchange of the meaning of matrix elements). 

2.3 Higher Hamiltonians 

Let us introduce "quasiclassic" hamiltonians: 



I k = TrA n L 1 (u)L 2 {u) . . . L k (u). 



Their classical limits in the sense of "Gr" are given by the following expressions: 
%\ = TrL(u)(n — 1)!; 

1 2 = (Tr 2 L(u) -TrL 2 (u))(n-2)\; (5) 

1 3 = (2TrL 3 (u) - 3TrL(u)TrL 2 (u) + Tr 3 L(u))(n - 3)! 



and can be taken as a basis of classical Gaudin hamiltonians. 
Consider formal expressions 



* = £(-i) j 'c?r, 

3=0 



*-3 



(6) 



where r = TrA n l = n\. 

Lemma 2 Sj /ias the following expansion on h 

Si {u) = h l QI t (u)+0(h l+1 ) 
where Qli is given by the formula 



QUiv) = Tr, n A n {L x {u) - d u )(L 2 (u) - d u ) . . . (L t (u) - d u )l 



(7) 



where Ai is the antisymmetrizer matrix in Mat^n an d Lj{u) is the quantum Lax operator 
(CJj lying in the j-th tensor component in Mat® l xn . 

Proof. Let us represent the formula for the generators of the Bethe subalgebra (jlj in 
a quite different manner: 



Ti(u - h)e~ thdu = TrA n e~" Uu T 1 {u)e~" Uu T 2 (u) ...<■' 



-hd u r 



-hd u 



TAu) 



where expressions on both sides are differential operators on u with values in U(qI u 
These formulas can be combined in a sort of a generating function: 

TrA^e-^T^u) - l)(e' hd -T 2 (u) - 1) . . . (e^T^u) - 1) 



u 



^ T .( M _ h )(_l)i-i C '| e 0'-0h^ 

3=0 



(9) 



Applying then both sides to the function 1 and using the fact that 

e- hd -T(u) - 1 = h(L(u) - d u ) + 0(h 2 ) 

we obtain 

Si {u) = \\ i TrA n {L x {u) - d u )(L 2 (u) - d u ) . . . {U{u) - d u )l + 0(W +1 ). 

□ 

Hence we have obtained that Qli commute. Indeed, [s»(it), Sj(v )] = due to the 
fact that Si(u) are linear combinations of generators of the Bethe subalgebra. Then, 
[QIi(u),QIj(v)] = in virtue of lemma H It is quite evident that the classical limit of 
Qli coincides with the classic hamiltonian, the element of the highest order in (J7J) is 

TrA n L 1 (u)L 2 (u) ...LAu). 



Thus we proved the main theorem of this paper: 
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Theorem 1 The quantities QIi(u) provide a quantization of the classical Gaudin model. 



Let us analyze the obtained quantum hamiltonians at small order: 

• 1 = 1 

s 1 = hTrL(u)(n- 1)! + 0(h 2 ); 
the first term at h is just the "semiclassic" hamiltonian I\(u). 

• 1 = 2 

s 2 = r 2 - 2n + to = h 2 (-cU>L( M )(n - 1)! + TrA n L 1 {u)L 2 {u)) + 0(h 3 ); 
the first term at h 2 can be expressed as 

h(u) - d u h(u). 

• 1 = 3 

s 3 = ^TrAn^Mu) - 9 tt Li(u)L2(u) - 2L 1 ( M )9 U L 2 ( M ) 
+ L 1 ( M )L 2 ( M )L 3 H) + 0(h 4 ). 

The first term at h 3 can be simplified as follows: 

h(u)-\duh{u) + d 2 Mu)\ 

where we have used the property TrA n d u L 1 (u)L 2 (u) = TrA n Li(u)d u L 2 (u) which is 
not true for the number of multipliers greater then 2. 

• 1 = 4 

s 4 = ^TrAni-dlL^u) + dlL^L^u) + 3Li(w)<9 2 L 2 (m) + ?>d u Li(u)d u L 2 (u) 

- <9 u Li(m)L 2 (m)L 3 (m) - 2L 1 (u)d u L 2 (u)L 3 (u) - 3Li(m)L 2 (m)<9 u L 3 (m) 

+ L 1 ( M )L 2 ( M )L 3 H)L 4 (M)) + 0(h 5 ). 

The first term here at h 4 also can be simplified: 

h(u) - d u h{u) + 2d 2 J 2 (u) - dlh(u) 
- TrA^dM^dMu) + L^duL^L^u) + 2L 1 (m)L 2 (m)9 u L 3 (m)) 

The arguments of the previous section give us that the first coefficients as like as their 
derivatives commute and we obtain that the quasiclassic hamiltonians for i = 1,...,3 
lie in our commutative family. However the higher hamiltonians i > 4 have unavoidable 
quantum corrections. 
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Concluding remarks 

Pull back to Y(gl n ) The proposed construction for quantum hamiltonians can be real- 
ized in a more general setup, it can be pulled back to the Yangian. In this case it provides 
a commutative family expressed only in the first order generators of the Yangian 

Qk(u) = T ri _ n A n (Tl l \u) - d u )(T^\u) -d u )... {Tf\u) - d u )l. 

In this way one could expect to quantize a quite wider class of models: rational matrices 
with higher order poles and the Hitchin system on singular curves. 

Bethe ansatz The obtained formulas for quantum hamiltonians provide a kind of de- 
composition which could be tractable in the context of the Bethe ansatz for higher di- 
mensional spin components of the Gaudin model. 
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